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INTRODUCTION  AND  BACKGROUND 


As  the  Air  Force  strives  to  maintain  its  excellence  in  science 
and  technology,  care  must  be  taken  to  investigate  unconventional 
areas  in  the  hopes  of  creating  a  scientific  breakthrough.  Such 
developments  keep  the  Air  Force  at  the  forefront  of  the  scientific 
community.  One  such  topic  is  the  zero-point  energy  (ZPE)  of  space. 
Theories  have  been  developed  regarding  how  to  tap  this  essentially 
infinite  supply  of  energy.  The  tapping  of  such  energy  could  have 
many  areas  of  application,  two  of  which  are  as  an  alternative 
energy  source  [Refs.  1,2]  and  as  a  field  generator  [Ref.  3].  The 
implication  of  nearly  limitless  energy  regardless  of  location  in 
space  points  to  an  obvious  utilization  of  the  ZPE  as  an  enabling 
means  of  previously  unrealized  types  of  propulsion. 

However,  research  and  development  regarding  the  ZPE  is  in  its 
infancy,  with  only  a  handful  of  researchers  having  done  serious 
work  in  the  field  [Ref.  3] .  In  this  paper,  a  theory  developed  by 
Nachamkin  will  be  investigated  [Ref.  4] .  His  theory  involves  the 
incidence  of  electromagnetic  (EM)  radiation  upon  a  dielectric 
sphere.  The  key  characteristics  of  this  phenomenon  are  a  sphere  of 
radius  a  and  propagation  constant  imbedded  in  an  infinite 
homogeneous  medium  of  constant  and  the  EM  wave  angular 
frequency,  co.  The  product  of  and  a  form  p,  the  parameter  of 
interest  for  this  type  of  interaction.  If  a  resonant  value  of  p 
can  be  identified,  then  it  would  be  desired  to  manufacture  two 
spheres  of  slightly  dissimilar  size,  but  with  each  able  to  attain 
the  same  value  of  p,  given  the  required  EM  radiation  angular 
frequency.  The  two  spheres  would  be  placed  within  close  proximity 
of  one  another  and  exposed  to  bombarding  EM  radiation  having  a 
range  of  frequencies  broad  enough  to  cover  the  difference  in 
resonant  O)  values  of  the  two  spheres.  In  turn,  the  resonant  waves 
established  would  work  in  conjunction  to  create  a  beat  frequency 
between  the  two  waves.  This  beat  frequency  could  be  rectified, 
thus  providing  a  source  of  electrical  energy  and  a  tap  of  the  ZPE. 

This  work  centered  upon  the  search  for  at  least  one 
identifiable  resonant  combination  of  ^sphere  radius  and  wave 
frequency,  as  mentioned  above.  Work  was  performed  at  the  Air 
Force's  Phillips  Laboratory  (PL)  at  Edwards  Air  Force  Base.  Use 
was  made  of  the  PL's  VAXCluster  for  computational  work. 
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THE  MIE  SCATTERING  AND  ITS  ASSUMPTIONS 


The  main  area  of  concern  addressed  in  this  report  is  the 
interaction  of  electromagnetic  radiation  with  a  dielectric  sphere; 
i.  e.,  the  diffraction  of  a  plane  wave  by  a  sphere,  more  commonly 
known  as  Mie  scattering  [Ref.  5]  .  It  is  assumed  that  the  sphere  is 
made  of  a  homogeneous  material  and  that  the  medium  surrounding  the 
sphere  is  a  vacuum.  The  incident  radiation  is  assumed  to  be  a 
plane  wave  propagating  in  the  z-direction.  Electrical  vibrations 
of  the  incident  wave  are  assumed  to  occur  in  the  x-direction,  with 
magnetic  vibrations  in  the  y-direction,  as  illustrated  in  Figure  1. 

As  explained  in  Stratton  [Ref.  6],  a  forced  oscillation  of 
free  and  bound  charges,  synchronous  with  the  applied  field,  arises 
when  a  periodic  wave  falls  incident  upon  a  body,  regardless  of  the 
sphere's  material.  This  creates  a  secondary  field  in  and  around 
the  body.  The  vector  sum  of  these  primary  and  secondary  fields 
gives  the  value  of  the  overall  field.  In  theory,  a  transient  term 
must  be  added  to  account  for  the  failure  of  the  boundary  conditions 
to  hold  during  the  onset  of  the  forced  oscillations.  However,  in 
practice  it  is  acceptable  to  consider  only  the  steady-state, 
synchronous  term  because  the  transient  oscillations  are  quickly 
damped  by  absorption  and  radiation  losses. 

A  plane  wave  falling  upon  a  sphere  is  the  simplest  of  such 
instances,  but  it  is  the  central  focus  of  this  study.  In 
Stratton's  notation,  the  sphere  has  a  radius  a  and  a  propagation 
constant  kj.  It  is  located  within  an  infinite,  homogeneous  medium 
with  propagation  constant  kj.  Region  1  and  region  2  have 
permittivities  Ei  and  Ej,  respectively,  and  permeabilities  and  jij, 
respectively.  Expanding  the  incident  field  in  vector  spherical 
wave  functions,  with  being  the  amplitude  of  the  time  harmonic 
electric  field,  yields  [Ref.  6]: 
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A  prime  represents  differentiation  with  respect  to  the  argument 
The  unit  vectors  ii,  £2,  and  13  represent  the  directions  of 
increasing  r,  9,  and  0,  respectively,  in  spherical  coordinates. 
They  are  presented  in  Figure  2 . 

The  electric  and  magnetic  fields  of  the  wave  transmitted  into 
the  sphere  (region  1),  i.  e.,  R  <  a,  can  be  expanded  using  the  same 
functions  [Ref.  6]: 
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If  j„(k3R)  is  replaced  by  h„'^‘ (kiR)  in  Eqs.  (3)  and  (4),  the 
functions  and  in  Eqs.  (5)  and  (6)  become  a‘^’  and  a‘^’.  The 
outgoing  fields  {R  >  a)  are  then  represented  by: 
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Figure  2 

Spherical  Coordinate  System  with  i  Unit  Vectors 
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It  should  be  noted  that  Eq.  (6)  is  replaced  with  k2  in  Eq. 
(8),  as  Eq.  (8)  represents  the  resultant  wave  in  the  medium 
surrounding  the  sphere  (region  2) . 

Being  sought  are  the  resonant  p  values  for  which  the  a,,'^  and 
b„*^  coefficients  are  infinite.  To  find  a^*^  and  bn*^,  the  boundary 
conditions  at  the  sphere  radius  (radius  R  =  a)  are  needed,  i.  e., 
continuity  of  the  E  and  H  values  at  the  surface.  Thus, 

X  (^j  +  E^)  =  i^x  (9) 


and 

X  (Hj  +  H^)  =  i^x  ,  (10) 


which  lead  to  two  pairs  of  inhomogeneous  equations: 
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Continuing  with  Stratton's  notation  [Ref.  6],  the  key  relations 
are : 


=  Nk^,  p  =  k^^a ,  k^a  =  Np  . 


(13) 


Spherical  Bessel  functions  of  the  first  kind  are  denoted  by  j„, 
while  those  of  the  third  kind  are  denoted  by  . 

Manipulation  of  Eqs .  (12)  and  (13)  gives  the  values  of  a,,-  and 
bn"  to  be: 
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At  a  resonance,  the  denominator  of  either  a„"  or  bn"  will  be  zero. 
A  FORTRAN  code  was  written  that  calculated  denominator  values  for 
varying  values  of  p  in  Eqs.  (14)  and  (15)  .  It  was  desired  that 
such  a  root  could  be  found  very  near  the  real  axis,  so  that  Real(p) 
»  Imaginary  (p)  .  Root  searches  were  started  by  using  the  code  to 
calculate  the  denominator  values  of  either  a„"  or  bn"  over  a  range 
of  p  values.  The  coefficient  type  ('a'  or  'b')  is  specified  in  the 
input  file. 

From  this  point  forward,  the  value  N  in  Eg.  (13)  will  be 
referred  to  as  KRATIO,  as  it  is  the  ratio  of  the  sphere  m.aterial's 
propagation  constant  kj  to  the  surrounding  medium's  propagation 
constant  k2.  This  is  done  in  order  to  avoid  confusion  with  the 
spherical  Bessel  fiinction  order  N  in  the  RES. FOR  FORTRAN  code, 
which  will  be  described  in  the  next  section. 
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The  program  is  utilized  to  find  p  values  that  correspond  to  a 
resonant  combination  of  angular  frequency  (0))  and  radius  (a)  for  a 
given  sphere  material  (region  1)  and  surrounding  medium  (region  2) . 
The  expression  for  p,  according  to  Stratton  [Ref.  6],  is: 


p  -  aicj  - 


(16) 


In  the  code,  p  is  identified  by  Fdi0_2,  as  it  corresponds  to  the 
surrounding  medium  (region  2)  .  RH0_1  (Pi)  corresponds  to  the 

sphere  material  (region  1)  and  is  given  by: 

Pi  -  (;ci/ic2)  p  .  (17) 


It  was  assumed  that  M’l  “  II2  =  M-o  =  47txl0"^  [H*m‘^]  and  that  62  =  Eq  = 
8.85419x10*^^  [F-m*^].  A  dielectric  constant  (Ej)  of  5.0  was  assumed 
for  the  sphere,  giving  the  relation:  E^  =  E^' E2  =  5.0- Ej.  Thus,  Ej 
was  considered  to  be  fully  real,  having  no  imaginary  part. 
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THE  CODE  (RES. FOR)  DESCRIPTION  AND  ITS  USAGE 


This  summary  describes  how  to  use  the  FORTRAN  code  RES. FOR. 
It  consists  of  the  Main  Program,  the  SUBROUTINE  SOLVE,  and  the 
SUBROUTINE  NEWTON.  Three  output  files  are  produced  by  RES.  They 
are  ITERRES.OUT,  REST. OUT,  and  RES2.0UT  .  One  input  file  is  used: 
RHORES.IN.  A  discussion  of  each  code  section  and  output  file 
appears  below,  followed  by  a  section  describing  how  to  use  the 
code.  Figure  3  is  a  flowchart  summarizing  the  code's  execution. 


MAIN  PROGRAM  AND  INPUT  FILE  (RHORES.IN) 

The  main  program  section  of  the  FORTRAN  code  RES. FOR  (Appendix 
A  contains  the  full  code)  begins  with  a  header,  followed  by  an 
OPTIONS  statement.  The  two  options  utilized  are  /G_FLOATING  and 
/I4.  Option  /G_FLOATING  was  utilized  in  the  program  to  allow  the 
handling  of  numbers  ranging  approximately  from  10'^°®  to  10*®°®  along 
both  the  real  and  imaginary  axes.  Option  /I4  defaults  all  integer 
variables  to  INTEGER*4  type.  All  variables  are  initially  declared 
as  COMPLEX*16,  followed  by  the  excepting  of  needed  INTEGER*4  and 
CHARACTER*!  variables  with  subsequent  declarations. 

Next  the  input  file,  RHORES.IN,  to  the  code  is  OPENed  and 
READ.  Examples  of  the  input  files  used  to  generate  the  needed 
results  are  found  in  Appendices  B,  C,  D,  and  E.  The  strategy  used 
in  varying  the  values  in  RHORES.IN  are  discussed  in  more  detail 
throughout  the  remainder  of  this  report.  DENTYP  is  the  type  of 
denominator  being  used  for  the  Newton-Raphson  iterations.  In 
Stratton  the  two  coefficients  are  an  and  bn,  so  here  DENTYP  is  a 
one-letter  variable  having  either  the  value  'a'  or  'b'  .  If  neither 
'b'  nor  'B'  is  placed  in  RHORES.IN,  the  default  is  'a'.  The 
default  IF-loop  follows  the  READing  of  RHORES.IN. 

NRITER  is  the  number  of  Newton-Raphson  iterations  to  be 
performed.  MU_1,  MU_2,  EPS_R,  and  EPS_2  are  as  explained  above. 
NBEG  and  NEND  give  the  range  of  spherical  'Bessel  function  orders  to 
be  calculated.  RHOMIN  and  RHOMAX  likewise  give  the  range  of  p 
values  to  be  tried  in  the  code,  with  RHODIV  being  the  number  of 
divisions  between  the  minimum  and  maximtim  values.  RHORES.IN  is 
then  CLOSEd,  and  ITERRES.OUT  is  OPENed.  All  of  the  input 
parameters  are  then  echoed  back  into  ITERRES.OUT.  A  message  tells 
the  user  the  name  of  the  three  output  files,  and  then  RESl.OUT  and 
RES2.0UT  are  OPENed.  KRATIO  is  then  calculated.  It  is  the  ratio 
of  to  kj,  as  seen  earlier. 

DO-loop  1100  steps  through  the  input  range  of  p  (RHO_2) 
values.  INIT  stores  the  value  of  RHO_2  before  Newton-Raphson 
recalculations.  RHO_l  is  found  as  mentioned  earlier,  and  ITER,  the 
current  Newton-Raphson  iteration,  is  set  to  0.  The  DOWHILE 
structure  is  used  for  the  iteration  loop,  with  the  condition  that 
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READ  in  RHORES.IN  Input  File 


'WRITE  Header  to/ 

iterres.out  / 

Output  File  / 


Begin  DO-loop  for  RHODIV  partitions  of  RHO_2  input  range 


STOP 


^iave  RHODIV  loop  pass^ 
^"^v^en  completed?^-''^ 


Calculate  RHO  1  and  RHO_2  values 


'WRITE  RHO_2  and  FRH( 
values  to  RES  1  .OUT 


Begin  DOWHILE  loop  for  number  of  iterations  required 


SUBROUTINE  SOLVE  calculates  Bessel 
function  values,  their  derivatives,  and  the 
a_n  and  b_n  denominator  values 


WRITE  RHO_2  and  appropriate 
denominator  value  to  / 

ITERRES.OUT  X 


WRITE  RHO_2,  function  orde^ 
and  inverse  denominator  / 
values  to  RES2,OUT  / 


SUBROUTINE  NEWTON  calculates 
next-guess  RHO_2  root  value 


Figurs  3 

Flowchart  for  RES. FOR  FORTRAN  Coda 
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ITER  must  be  less  than  or  equal  to  NRITER.  The  first  step  in  tihe 
DOWHILE  loop  is  to  call  the  SUBROUTINE  SOLVE  in  order  to  determine 
the  denominator  values  of  a^  (DENA)  and  (DENB) ,  as  well  as  the 
derivatives  of  DENA  and  DENB,  DENA$  and  DENB$,  respectively.  it 
should  be  noted  that  '$'  represents  a  first  derivative  of  the 
variable  with  respect  to  p,  while  '$$'  represents  a  second 
derivative.  Depending  on  the  type  of  denominator  (aj,  or  b^) 
currently  being  scrutinized,  RH0_2  and  the  appropriate  DEN  term  are 
written  to  ITERRES.OUT.  Newton-Raphson  is  then  called  upon  to 
determine  the  new  RH0_2  value  for  the  next  iteration,  passing  the 
appropriate  DEN  and  DEN$  values.  RHO_l  is  then  recalculated,  and 
ITER  is  stepped  one. 

The  beginning  RH0_2  value  and  FRHO  (a  function  of  RHO_2)  are 
then  written  to  RESl.OUT  for  the  current  p  division  of  the  1100  DO- 
loop.  Finally,  the  three  output  files  are  CLOSEd  and  code 
operations  are  stopped. 


SUBROUTINE  SOLVE 

This  SUBROUTINE  represents  the  largest  portion  of  the  code. 
In  it,  the  Jn  and  spherical  Bessel  function  values  up  to  the 
needed  order  N  are  calculated.  These,  in  turn,  yield  the  values 
for  DENA,  DENB,  DENA$ ,  DENB$ ,  and  FRHO.  Four  1-dimensional  arrays 
are  used  in  the  calculation  of  the  j's  and  y's.  By  having 
capacities  of  10000,  these  arrays  may  contain  j  and  y  values  up  to 
the  9995th  order.  It  is  not  a  full  10000  because  array  indexes  of 
1,  2,  9999,  and  10000  are  required  in  the  calculations  of  second 
derivative  ($$)  terms,  and  an  N=0  index  must  be  included.  The 
indexes  (shown  in  parentheses  below)  are  offset  by  3  from  the  true 
function  order  (N) : 


(1)  — -> 

N=-2 

(2)  ---> 

N=-l 

(3)  ---> 

N=  0 

(4)  — -> 

N=+l 

(5)  ~  — > 

N=+2 

(9998)  — -> 

N=+9995 

The  maximum  order  index  (MAXORD  and  MAXOR2)  must  be  offset  2 
because  of  the  needed  allocation  space  for  the  N=-2  and  N=-l  cases. 
MAXORD  will  never  have  a  value  less  than  10,  as  given  by  the  JMAXO 
functions.  MAXORD  is  the  upper  order  limit  for  RHO_l-only 
functions,  while  MAXOR2  is  the  same  for  fxinctions  of  RHO_2  only. 
The  Jn's  for  RHO_l  (found  in  the  array  JNRHOl)  and  RHO_2  (array 
JNRH02)  are  then  solved  as  was  done  by  Nachamkin  in  his  program 
TESTJY.FOR.  This  is  followed  by  the  calculations  of  the  y„'3  for 
RHO_2  (array  YNRH02) .  No  values  are  needed  for  RHO_l  because 
the  terms  in  the  denominators  are  j„'s  (spherical  Bessel  functions 
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of  the  1st  kind),  which  are  functions  of  RH0_1  only.  The 
(3rd  kind)  are  functions  of  RH0_2  only  and  are  equivalent  to: 


(18) 


Upon  completion  of  the  y„  calculations,  some  of  the  variables 
are  initialized  and  the  current  RH0_2  value  is  written  to  RES2.0UT. 
DO-loop  4000  then  finds  all  of  the  denominator  and  derivative 
values  based  on  recurrence  relations  found  in  a  handbook  of 
mathematical  functions  [Ref.  7].  A,  B,  C,  and  CINV  are  used  to 
eliminate  some  of  the  bulk  in  the  recurrence  relations  and  make 
them  more  readable  and  more  easily  editted.  The  4000  DO-loop  uses 
J  as  its  counting  variable,  starting  at  0  and  going  to  MAXN,  as 
found  in  the  y„  calculations.  ZERO  represents  the  index 
corresponding  to  the  actual  Jth  function  being  observed. 

ON  and  DN  are  the  inverses  of  the  an  and  bn  denominator  values, 
respectively.  SQCN  and  SQDN  are  simply  the  squares  of  the  absolute 
values  of  CN  and  DN,  respectively.  FRHO  represents  the  running  sum 
of  the  SQCN  and  SQDN  terms  combined  as  the  loop  is  passed  through 
MAXN+1  times.  SQCNTOT  and  SQDNTOT  provide  a  more  detailed  account 
of  the  gradual  accxomulation  of  the  FRHO  term  through  the  separate 
running  sxjiranations  of  the  SQCN  and  SQDN  values.  A  special  note  is 
that  the  second  derivatives  of  j„  and  have  been  commented  out 
because  the  Newton-Raphson  method  seemed  to  converge  better  in  some 
cases  without  these  terms  included.  This  should  be  investigated 
further,  however,  before  leaving  these  terms  out  completely.  The 
second  derivative  terms  are  only  utilized  in  the  Newton-Raphson 
portion  of  the  code,  so  leaving  them  out  is  not  detrimental  to  the 
performance  of  the  other  portions  of  the  code.  This  is  not 
necessary  a  limitation  of  the  code,  but  rather  a  possible 
limitation  of  the  use  of  the  Newton-Raphson  method  for  functions 
involving  spherical  Bessel  functions. 

If  the  number  of  Newton-Raphson  iterations  wanted  is  0  and 
RHO_2  is  being  held  constant,  then  the  current  J  value  and  all 
SQCN,  SQDN,  SQCNTOT,  SQDNTOT  values  will  be  written  to  RES2.0UT  . 
This  will  also  occur  if,  in  addition  to  0  Newton-Raphson  iterations 
being  input,  a  single  order  is  being  scrutinized  for  a  resonance 
(NBEG=NEND)  and  J  has  reached  that  order's  value  ( J=NBEG=NEND) . 


SUBROUTINE  NEWTON 

In  this  prograun  structure,  the  method  of  Newton-Raphson  is 
used  to  attempt  to  converge  on  a  root  in  either  the  a„  or  b^ 
denominator  equations,  depending  on  which  DENTYP  is  passed  to  it. 
A  precaution  has  been  taken  to  let  the  user  know  in  the  unlikely 
case  that  DEN$  equals  (0.D0,0.D0),  which  immediately  would  kill 
program  execution  due  to  division  by  zero.  After  determining  the 
new  RHO_2  value,  control  returns  to  the  Main  Program. 
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OUTPUT  FILE  #1  (ITERRES.OUT) 

See  Appendix  E. 


OUTPUT  FILE  #2  (RESl.OUT) 

See  Appendices  B  and  D. 

OUTPUT  FILE  #3  (RES2.0UT) 

See  Appendix  C. 


RES. FOR  CODE  USAGE 

Earlier  work  by  Nachamkin  [Ref.  4]  had  yielded  possible 
neighborhoods  of  root  locations  along  the  Real  axis.  Thus,  these 
prospective  resonance  locations  were  checked  using  the  current 
code.  It  was  found  that  the  current  code's  findings  supported 
Nachamkin' s  earlier  results.  The  checking  of  the  possible  root 
locations  began  by  inputting  a  range  of  p  values  in  the 
neighborhood  of  the  prospective  root.  The  range  was  divided  into 
a  specified  number  of  intervals.  For  each  such  interval,  the 
inverse  value  of  the  denominator  was  determined. 

The  code  is  utilized  in  four  ways  during  the  determination  of 
a  root.  First,  a  range  of  RH0_2  values  is  studied  to  find  any  FRHO 
peaks  in  that  range.  Second,  once  a  peak  has  been  chosen,  the 
function  order  n  giving  the  dominant  FRHO  term  is  determined.  This 
also  gives  a  clue  as  to  whether  the  peak  is  due  to  a  magnetic 
resonance  (an  approaches  infinity)  or  an  electrical  resonance  (bn 
approaches  infinity)  .  If  the  dominant  FRHO  value  comes  from  tbe  CN 
term,  one  uses  'a'  (for  DENTYP)  in  RHORES.IN,  while  'b'  is  used  for 
a  dominant  DN  term.  Third,  a  large  number  of  Newton-Raphson 
iterations  is  performed  in  order  to  converge  upon  a  root  RH0_2 
value.  A  hundred  or  more  iterations  is  possible.  This  root  may  or 
may  not  lie  near  the  initial  RHO_2  value  wWich  gave  the  peak  in  the 
FRHO  function.  Finally,  after  finding  a  root,  the  root  is  checked 
to  make  sure  it  has  a  dominant  term  for  the  same  order  n  originally 
determined.  Then  the  root  is  used  in  the  code  as  the  initial  RH0_2 
value  to  see  if  the  denominator  does,  indeed,  converge  toward  0.0. 

In  order  to  understand  the  process  better,  the  root  identified 
in  this  study  will  be  used  as  an  example.  The  spherical  Bessel 
function  order  N  varies  from  0  to  100,  because  the  resonant  order 
has  not  yet  been  determined.  The  real  portion  of  RHO_2  is  held 
constant,  while  the  complex  portion  varies  from  -1.0  to  0.0  with 
100  divisions.  In  Appendix  B,  the  RESl.OUT  file,  one  sees  that  the 
FRHO  value  achieves  its  maximum  value  for  the  -0.63  imaginary  mark. 
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since  the  peak  value  is  now  determined  to  two  significant 
figures,  one  may  determine  the  N  value  that  contributes  the  most  to 
FRHO,  i.  e.,  the  resonant  spherical  Bessel  function  order.  This  is 
seen  in  Appendix  C,  which  displays  the  contents  of  the  output  file 
used  for  this  purpose--RES2 .OUT.  The  RH0_2  value  is  held  constant, 
while  the  value  of  N  varies  from  0  to  100.  For  N  =  58,  the  DN 
value  (resulting  from  b^'s  denominator)  is  largest  by  about  four 
orders  of  magnitude  over  all  other  N  values.  Thus,  one  sees  that 
the  N  range  may  now  be  narrowed  down  to  a  single  order,  58,  for 
determination  of  the  root.  Also,  DENTYP  would  be  'b'  because  the 
spike  in  the  FRHO  value  results  from  the  bn  denominator. 

Having  narrowed  down  the  n  value  to  a  single  order,  further 
RHO_2  division  runs  are  made.  The  final  run  is  shown  in  Appendix 
D,  which  is  similar  to  Appendix  B,  but  now  is  based  on  one  N  value 
and  is  much  more  accurate.  This  time  the  FRHO  value  reaches  order 
of  magnitude  10^®,  resulting  from  a  complex  portion  of  p  equal  to 
-0.634786707197.  The  full  1000  RHO_2  divisions  have  been  narrowed 
down  to  approximately  250  divisions,  located  in  Appendix  D  as  well. 

Finally,  this  root  value  may  be  checked  using  the  RHORES.IN 
input  file  in  Appendix  E.  Twenty  iterations  are  calculated,  with 
N  and  RHO_2  held  constant.  To  observe  the  iterations  of  the  NEWTON 
SUBROUTINE,  the  ITERRES.OUT  output  file  is  used.  Appendix  E  also 
contains  the  resulting  ITERRES.OUT  file. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


The  one  major  root  of  p  which  was  identified  had  a  value  of: 

REAL(p)  =  +66.39752607619131 
IMAGINARY(p)  =  -0.6347867071968998  . 

Thus,  one  may  determine  a  possible  aco  combination  which  would  have 
this  root  value.  For  p,  e  =  Eq  and  |i  =  [Xq, 

p  -  -  aa>/c  .  (19) 


In  SI  units,  the  speed  of  light  c  =  2.99792458x10®  [m/s].  If  an  a 
value  of  10‘®  [m]  is  assumed,  then: 


(0  -  -££  -  l.9919xl0i«  -  i  1.9044x10^*  [rad/a]  .  (20) 


This  then,  is  the  angular  frequency  required  within  the  impingent 
EM  radiation  in  order  to  create  a  resonant  situation.  Other 
resonances  were  indicated  by  the  outputs  from  the  code,  but  this 
was  the  most  prevalent  one. 

Of  course,  no  complex-frequency  plane  waves  exist.  Thus,  the 
code  was  used  by  considering  only  the  real  portion  of  the  above 
root,  setting  the  imaginary  portion  equal  to  zero.  Upon  doing  this, 
however,  the  code  became  insensitive  to  any  root  in  the  vicinity  of 
the  root's  real  portion.  Therefore,  a  number  of  runs  looking  at 
ranges  of  real-only  p  values  in  the  p  =  1  to  100  neighborhood  were 
executed.  The  results  of  the  p  =  40  to  50  range  are  found  in  Figure 
4,  showing  FRHO  values  versus  the  real  portion  of  p.  Possible 
resonances  are  indicated  by  sharp  rises  in  the  Figure  4  plot.  One 
such  resonance  is  detailed  more  closely  in  the  Figure  5  graph, 
showing  an  approximate  two  order  of  magnitude  rise  over  the 
surrounding  curve  region  at  a  RHO_2  (p)  value  of  approximately 
+43.51883  .  Yet  many  elusive  resonances  surely  remain,  awaiting 
more  code  executions  targeted  along  the  real  axis. 

Much  work  still  remains  in  finding  more  resonances  and  in 
studying  other  areas  of  the  theory.  A  source  of  EM  radiation  having 
a  broad  enough  range  of  frequencies  to  achieve  resonances  between 
two  chosen  spheres  needs  to  be  selected.  Then,  one  should  analyze 
the  beat  frequency  produced  by  the  interaction  of  the  two  resonant 
waves,  as  well  as  the  effect  of  separation  distance  of  the  two 
spheres  on  the  beat  frequency.  Finally,  a  method  of  rectifying 
this  beat  frequency  should  be  established  using  currently  available 
equipment,  if  possible.  It  is  also  important  to  know  how  much 
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energy  is  available  at  the  resonant  points.  As  a  practical  matter, 
manufacturing  processes  must  be  investigated  that  would  alio 
structures  to  be  fabricated  with  close  enough  tolerances  to  be  o 
use . 


FRHO 


FRHO  vs.  REAL(RHO_2) 


Figure  4 

FRHO  vs.  REAL(RHO_2) 
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FRHO 


FRHO  vs.  REAL(Resonant  RHO_2) 


REAL(RHO_2) 


Figure  5 

FRHO  Around  RHO_2  Resonant  Value 
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APPENDIX  A 


RES. FOR 


PROGRAMMER:  LARRY  T.  COX,  JR. 

ORGANIZATION:  OLAC-PL/RKFE 

ADDRESS:  EDWARDS  AFB,  CA  93523-5000 
DATE  OF  ORIGIN:  4  AUGUST  1992 
LATEST  REVISION:  6  AUGUST  1992 

OPTIONS  /G_FLOATING  /I4 

C23456 

PROGRAM  RESONATOR 
IMPLICIT  COMPLEX*16  (A-Z) 

INTEGER*4  I,  ITER,  J, K,  L, M, N,  WRITER, RHODIV, NBEG, NEND, MAXN 
CHARACTER*!  DENTYP 

OPEN  (UNIT=1, FILE='RHORES. IN' ,FORM=' FORMATTED' , STATUS='OLD' , 

&  ACCESS=' SEQUENTIAL' ) 

READ  (1,99)  DENTYP 
99  FORMAT  (Al) 

READ  (1,*)  WRITER  !  Number  of  Newton-Raphson  iterations 

READ  ( 1 , * )  MU_1 

READ  (1,*)  EPS_R 

READ  (1 ,  * )  MU_2 

READ  (1,*)  EPS_2 

READ  (1,*)  NBEG  !  Function  of  the  Nth  kind 

READ  (1,*)  NEND  !  Function  of  the  Nth  kind 

READ  (1,*)  RHOMIN  !  Initial  value  of  RHO_2 

READ  ( 1 ,  * )  RHOMAX 
READ  (1,*)  RHODIV 
CLOSE  (UNIT=1) 

EPS_1  =  EPS_R  *  EPS_2 

IF  (DENTYP. EQ.'B' .OR. DENTYP. EQ.'b' ^  THEN 
DENTYP  =  'b' 

ELSE 

DENTYP  =  'a' 

ENDIF 

OPEN  (UNIT=1,FILE=' ITERRES .OUT' , FORM=' FORMATTED' , STATUS='NEW' , 
&  ACCESS®' SEQUENTIAL' ) 

WRITE  (1,*)  DENTYP, '_n  denominator  roots' 

WRITE  (1,*)  'mu_l  =  ',MU_1 
WRITE  (1,*)  'mu_2  =  ',MU_2 
WRITE  (1,*)  'epsilon_l  =  ',EPS_1 
WRITE  (1,*)  'epsilon_2  =  ',EPS_2 

WRITE  (1,*)  'Functions  of  the  ',NBEG, 'th  to  ',NEND,'th  kind.' 
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1095 


1098 

1099 


WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 


PRINT 

PRINT 

PRINT 


(1,  *) 

(1,*)  '  RE(RHO)  IM(RHO)' 

(1,*)  DREAL (RHOMIN) , DIMAG (RHOMIN) 

(1,*)  DREAL (RHOMAX) , DIMAG (RHOMAX) 

(1,*)  RHODIV, '  divisions' 

(1,*) 

(1,*)  '  Re(rho)  Im(rho)', 

'  Re (Denom  ' , DENTYP, '_n)  Im(Denom 

DENTYP, '_n) ' 

'Output  file  containing  iterations  is  ITERRES.OUT 
*,  'Other  output  files  are:' 

*,  'RESl.OUT  and  RES2.0UT  .' 


/ 


OPEN  (UNIT=2,FILE=' RESl.OUT' , FORM=' FORMATTED' , STATUS='NEW' , 
ACCESS=' SEQUENTIAL' ) 

OPEN  (UNIT=3, FILE='RES2. OUT' ,FORM=' FORMATTED' , STATUS='NEW' , 
ACCESS=' SEQUENTIAL' ) 

KRATIO  =  CDSQRT (MU_1*EPS_1  /  (MU_2*EPS_2) ) 

WRITE  (2,*)  'Rho''s  with  Functions  of  Rho''s' 

DO  1100  L  =  l,RHODIV 

RHO_2  =  RHOMIN  +  DCMPLX (L-l) * (RHOMAX-RHOMIN) /DCMPLX (RHODIV) 

INIT  =  RHO_2 

RHO_l  =  KRATIO  *  RHO_2 

ITER  =  0 

DOWHILE  (ITER.LE.NRITER) 

CALL  SOLVE (KRATIO, RHO_l , RHO_2 , DENA, DENA$, N, RHODIV, 

MU_1 , MU_2 , DENB, DENB$ , NBEG, NEND , FRHO, NRITER 

IF  (DENTYP .EQ. 'b' )  THEN 

WRITE  (1,1095)  RHO_2,DENB 

ELSE 

WRITE  (1,1095)  RHO_2,DENA 

END  IF 

FORMAT  (2E24 .16,2E16.8) 

IF  (DENTYP  .EQ. 'b' )  THEN' 

CALL  NEWTON (ITER, NRITER, DENB, DENB$, RHO_2) 

ELSE 

CALL  NEWTON (ITER, NRITER, DENA, DENA$ , RHO_2 ) 

END  IF 

RHO_l  =  KRATIO  *  RHO_2 
ITER  =  ITER  +  1 

ENDDO 

FORMAT  {5E16.8) 

WRITE  (2,1099)  INIT, FRHO  , 

FORMAT  (E24.16,E24.16,E24.16,E8.2) 

WRITE  (1,*) 
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o  o  o 


1100  CONTINUE 


CLOSE  (UNIT=1) 
CLOSE  (UNIT=2) 
CLOSE  (UNIT=3) 

STOP 

END 


C23456 

OPTIONS  /G_FLOATING  /I4 

SUBROUTINE  SOLVE (KRATIO, RHO_l , RHO_2 , DENA, DENA$, N, RHODIV, 

&  MU_1 , MU_2 , DENB , DENB$ , NBEG , NEND , FRHO , NRI TER) 

IMPLICIT  COMPLEX*16  (A-Z) 

INTEGER*4  I, lORD, J,  K, L, M, MAXORD, MAXOR2, N, ORDER, RHODIV, ZERO 
INTEGER* 4  NBEG, NEND, MAXN, NRI TER 
REAL*8  SQCN, SQDN, A, B, C, CINV, SQCNTOT, SQDNTOT 
COMPLEX*16  JNRHOl (10000) , JNRH02  (10000)  ,  YNRH02 (10000) 

COMPLEX* 1 6  H1NRH02 (10000) 

DATA  JNRHOl  /lOOOO  *  (0.D0,0.D0)/ 

DATA  JNRH02  /lOOOO  *  (0.D0,0.D0)/ 

DATA  yNRH02  /lOOOO  *  (0.D0,0.D0)/ 

DATA  H1NRH02  /lOOOO  *  (0.D0,0.D0)/ 

ORDER  =  NEND+1 

MAXORD  =  2  +  JIDINT(DREAL(RHO_l) )  + 

&  JIDINT(DREAL(CDSQRT( ( 150 . DO , 0 . DO ) *RHO_l) ) ) 

MAXORD  =  UMAX 0 (MAXORD, 10) 

MAXOR2  =  2  +  JIDINT (DREAL (RHO_2) )  + 

&  JIDINT (DREAL (CDSQRT ( (150 .DO, 0 .DO) *RHO_2) ) ) 

MAXOR2  =  JMAXO (MAXOR2, 10) 

IF  (MAXORD. GT. 9999)  THEN 

PRINT  *,  '  RHO_l  TOO  BIG  !!!''(=  DREAL  (RHO_l)  ,')  ' 

STOP 

ELSEIF  (MAXOR2.GT.9999)  THEN 

PRINT  *,  '  RHO_2  TOO  BIG  !!!  (=  DREAL (RHO_2) ,') ' 

STOP 

END  IF 

Begin  the  jn' s  for  RHO_l 

RHOIINV  =  (l.D0,0.D0) /RHO_l 
RH012IN  =  RHOIINV  +  RHOIINV 

B1  =  -RHOIINV  -  RH012IN 

JNRHOl (1)  =  -CDCOS (RHO_l) /RHO_l**2  -  CDSIN (RHO_l) /RHO_l 
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o  o  o 


JNRH01(2)  =  CDCOS (RH0_1) /RH0_1 
RIO  =  -CDSIN (RHO_l) /RHO_l 
Rll  =  (RIO+CDCOS (RHO_l) ) *RH01INV 
BETAl  =  B1 

JNRHOlO)  =  RlO/BETAl 
lORD  =  4 

B1  =  B1  -  RH012IN 

YNRH02 (lORD)  =  ( 1 . DO , 0 . DO ) /BETAl 

BETAl  =  B1 

JNRHOl(IORD)  =  Rll/BETAl 

DO  2000  lORD  =  5,MAXORD+4 
B1  =  B1  -  RH012IN 
YNRH02(I0RD)  =  (1 .DO, 0 .DO) /BETAl 
BETAl  =  B1  -  YNRH02(I0RD) 

JNRHOl(IORD)  =  -JNRHOl (IORD-1) /BETAl 
2000  CONTINUE 

DO  2100  lORD  =  MAXORD+3, 3, -1 

JNRHOl (lORD)  =  JNRHOl (lORD) -YNRH02 (IORD+1) *JNRH01 (IORD+1) 
2100  CONTINUE 

DO  2200  lORD  =  MAXORD+3, 3, -1 

JNRHOl (IORD+1)  =  JNRHOl (lORD) 

2200  CONTINUE 

JNRHOl (3)  =  -RIO 

Begin  the  jn' s  for  RHO_2 

RH02INV  =  (1 .DO, 0 .DO) /RHO_2 
RH022IN  =  RH02INV  +  RH02INV 

B2  =  -RH02INV  -  RH022IN 

JNRH02(1)  =  -CDCOS (RHO_2) /RHO_2**2  -  CDSIN (RHO_2) /RHO_2 

JNRH02(2)  =  CDCOS (RHO_2) /RHO_2 

R20  =  -CDSIN (RHO_2) /RHO_2 

R21  =  (R20+CDCOS (RHO_2) ) *RH02INV  " 

BETA2  =  B2 

JNRH02(3)  =  R20/BETA2 
lORD  =  4 

B2  =  B2  -  RH022IN 

YNRH02(IORD)  =  (1 .DO, 0 .DO) /BETA2 

BETA2  =  B2 

JNRH02(I0RD)  =  R21/BETA2 

DO  2300  lORD  =  5,MAXOR2+4 
B2  =  B2  -  RH022IN 
YNRH02(I0RD)  =  (1 .DO, 0 .DO) /BETA2 
BETA2  =  B2  -  YNRH02(I0RD) 
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JNRH02(I0RD)  =  - JNRH02 ( IORD-1 ) /BETA2 
2300  CONTINUE 

DO  2400  lORD  =  MAXOR2+3, 3, -1 

JNRH02(I0RD)  =  JNRH02 (lORD) -YNRH02 (IORD+1) *JNRH02 (IORD+1) 
2400  CONTINUE 

DO  2500  lORD  =  MAXOR2+3, 3, -1 

JNRH02 (IORD+1)  =  JNRH02 (lORD) 

2500  CONTINUE 

JNRH02(3)  =  -R20 

Begin  the  yn's  for  RHO_2 

yNRH02 ( 1 )  =  - JNRH02 ( 4 ) 

YNRH02(2)  =  JNRH02(3) 

YNRH02(3)  =  -CDCOS (RHO_2)  /  RHO_2 

YNRH02(4)  =  -JNRH02(3)  +  YNRH02 (3)  /  RHO_2 

lORD  =  4 

TN2  =  RH02INV 

TN2P  =  (l.D0,0.D0) 

IF  (NEND.GT.MAXORD)  THEN 
MAXN  =  MAXORD 

ELSE 

MAXN  =  NEND 

END  IF 

DOWHILE  ( lORD . LE . (MAXORD+4 ) ) 
lORD  =  IORD+1 
TN2  =  TN2  +  RH022IN 

YNRH02(IORD)  =  TN2*YNRH02 ( IORD-1)  -  YNRH02 (IORD-2) 

TN2P  =  TN2  +  RH022IN 

ENDDO 

DO  3000  I  =  1, MAXORD+4 

H1NRH02(I)  =  JNRH02(I)  +  (0 .DO, 1 .DO) *YNRH02 (I) 

3000  CONTINUE 

ONE  =  (l.D0,0.D0) 

TWO  =  (2. DO, 0. DO) 

FRHO  =  (0.D0,0.D0) 

SQCNTOT  =  O.DO 
SQDNTOT  =  O.DO 
WRITE  (3,*)  'RHO  =  ',RHO_2 

DO  4000  J  =  0,MAXN 
A  =  DFLOAT(J) 

B  =  A  +  l.DO 
C  =  2.D0*A  +  l.DO 
CINV  =  l.DO  /  C 
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ZERO  =  J+3 

JN  =  JNRHOl(ZERO) 

JN$  =  {A*JNRH01(ZER0-1)  -  B* JNRHOl (ZERO+1 ) )  *  CINV 

JNM1$  =  (A* JNRHOl {ZERO-2)  -  B*JN)  *  CINV 
JNP1$  =  (A*JN  -  B*JNRH01(ZER0+2) )  *  CINV 
HIN  =  H1NRH02 (ZERO) 

H1N$  =  (A*H1NRH02 (ZERO-1)  -  B*H1NRH02 (ZERO+1 ) )  *  CINV 
H1NM1$  =  (A*H1NRH02 (ZERO-2)  -  B*H1N)  *  CINV 
H1NP1$  =  (A*H1N  -  B*H1NRH02 (ZERO+2) )  *  CINV 

C  JN$$  =  (A*JNM1$  -  B*JNP1$)  *  CINV 

C  H1N$$  =  (A*H1NM1$  -  B*H1NP1$)  *  CINV 

DENA  =  -  HIN  *  (RHO_l*JN$+JN) 

+  JN  *  (RH0_2*H1N$+H1N) 

DENA$  =  +  JN*  (RHO_2*HlN$$  +  2  .D0*H1N$)  +  .JN$*  (RH0_2*H1N$+H1N) 

-  HIN* (RHO_l*JN$$+2 .D0*JN$) 

-  H1N$* (RHO_l*JN$+JN) 

DENB  =  +  HIN  *  (RHO_l*JN$+JN) 

-  KRATIO**2*JN  *  (RH0_2*H1N$+H1N) 

DENB$  =  +  HIN* (RHO_l*JN$$+2 .D0*JN$) 

+  H1N$* (RHO_l*JN$+JN) 

-  KRATIO**2*JN* (RHO_2*HlN$$+2 .D0*H1N$) 

-  KRATIO**2*JN$* (RH0_2*H1N$+H1N) 

CN  =  ONE  /  DENA 
DN  =  ONE  /  DENB 

SQCN  =  (DREAL (CN) ) **2  +  (DIMAG (CN) ) **2 
SQDN  =  (DREAL(DN) ) **2  +  (DIMAG (DN) ) * *2 
FRHO  =  FRHO  +  SQCN  +  SQDN 
SQCNTOT  =  SQCNTOT  +  SQCN 
SQDNTOT  =  SQDNTOT  +  SQDN 

IF  (NRITER.EQ.O.AND.RHODIV.EQ.l)  THEN 
WRITE  (3,*)  'ORDER  =  ',J 
WRITE  (3,3995)  SQCN, SQCNTOT, SQDN, SQDNTOT 
ELSEIF  (NBEG.EQ.NEND.AND. J.EQ.NEND.AND.NRITER.EQ.O)  THEN 
WRITE  (3,*)  'ORDER  =  ',J 
WRITE  (3,3995)  SQCN, SQCNTOT, SQDN, SQDNTOT 
WRITE  (3,*) 

ENDIF 

3995  FORMAT  (4G20.10) 

4000  CONTINUE 

RETURN 

END 


C23456 

OPTIONS  /G_FLOATING  /I4 

SUBROUTINE  NEWTON (ITER, NRITER, DEN, DEN$, RHO_2) 
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IMPLICIT  C0MPLEX*16  (A-Z) 


5097 


INTEGER*4  ITER,NRITER 

IF  (DEN$.EQ.  (0.D0,0.D0)  )  THEN 
ITER  =  NRITER 

WRITE  (1,5097)  ' DEN$=  (0 .DO, 0  .DO) ' 
FORMAT  (A) 

ELSE 

DELRHO  =  -DEN  /  DEN$ 

RHO_2  =  RHO_2  +  DELRHO 

ENDIF 

RETURN 

END 
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APPENDIX  B 


RHORES.IN  and  RES 1. OUT  (1st  run) 


b 

0 

(1 .25663706144D-6, 0 .DO) 
(5 .DO, 0 .DO) 

(1 .25663706144D-6, 0 .DO) 
(8 . 85418782D-12, 0 .DO) 


0 

100 

(66.3975260761913D0, -1 .DO) 
(66.3975260761913D0, O.DO) 
100 
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Rho ' 3  with  Functions 
0 . 663975:507619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0 . 6639752607619130E+02  - 
0.6639752607619130E+02  - 
0 . 6639752607619130E+02  - 
0 . 6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0 . 6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0. 6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02 
0.6639752607619130E+02 
0.6639752607619130E+02 
0.6639752607619130E+02 
0.6639752607619130E+02 
0.6639752607619130E+02 
0.6639752607619130E+02 
0.6639752607619130E+02 
0.6639752607619130E+02 
0.6639752607619130E+02 
0.6639752607619130E+02 
0.6639752607619130E+02 
0.6639752607619130E+02 
0.6639752607619130E+02 
0.6639752607619130E+02 
0.6639752607619130E+02 
0.6639752607619130E+02 
0.6639752607619130E+02 
0.6639752607619130E+02 
0.6639752607619130E+02 
0.6639752607619130E+02 


;  o  f  Rho ' s 

0 . lGCOOOOCE+01 

0.9900000CE+00 

0 . 98000000E+00 

0.9700000CS+00 

0.9600000CE+00 

0.9500000CE+00 

0.94000000E+00 

0.9300000CS4.00 

0.92000000E+00 

0.9100000CE+00 

0.9000000GE+00 

0.89000000E+00 

0.8800000CE+00 

0.87000000E+00 

0.8600000CE+00 

0. 85000000E+00 

0.8400000CE+00 

0.83000000E+00 

0.8200000CE+00 

0.8100000CE+00 

0.80000000E+00 

0.7900000CE+00 

0.78000000E+00 

0.7700000CE+00 

0.7600000CS+00 

0.75000000E+00 

0.7400000CE+00 

0.7300000CE+00 

0.72000000E+00 

0.71000000E+00 

0.70000000E+00 

0.69000000E+00 

0.68000000E+00 

0.67000000E+00 

0.6600000CE+00 

0.65000000E+00 

0.64000000E+00 

0.63000000E+00 

0.62000000E+00 

0.61000000E+00 

0.60000000E+00 

0.59000000E+00 

'0.58000000E+00 

■0.57000000E+00 

■0.56000000E+00 

•0.55000000E+00 

•0.54000000E+00 

■0.53000000E+00 

•0.52000000E+00 

■0.51000000E+00 

•0.50000000E+00 

•0.49000000E+00 

■0.48000000E+00 

■0.47000000E+00 

■0.46000000E+00 

•0.45000000E+00 

■0.44000000E+00 

-0.4300000CE+00 

■0.42000000E+00 

■0.41000000E+00 

■0.40000000E+00 

■0.39000000E+00 

-0.38000000E+00 


0.583294601 9812745 E+C5 

0.5942793342897465E+C5 

0 . 6065642728623391E^05 

0 . 6202102098364594E+05 

0.6352919605634567E^-05 

0.6518974139281166S+05 

0.6701302511900682E+05 

0.6901125835360250S+05 

0.7119881849116909E+05 

0.7359264829137213E+05 

0.7621275203192639E+05 

0.7908281686893944E+05 

0 .8223099716105772E+05 

0.8569091309061948E+05 

0.8950293435562658E+05 

0 . 9371584798417194E+05 

0.9838905117578449E+05 

0.1035954732316571E+06 

0.1094255279987954E+06 

0.1159925519140877E+06 

0.1234404315763505E+06 

0.1319545392892882E+06 

0. 1417778071454112E+06 

0.1532350356617195E+06 

0.1667708665727319E+06 

0.1830113327824183E+06 

0.2028679161262282E+06 

0.2277221518793277E+06 

0.2597718390884264E+06 

0.3027240302222338E+06 

0.3632952876402438E+06 

0.4547938649770824E+06 

0.6068445954717640E+06 

0.8969560608025888E+06 

0.1584843084218144E+07 

0.4032547676841964E+07 

0.3299221352934197E+08 

0.3945092293223719E+08 

0.4389945914409366E+07 

0.1737334773076357E+07 

0.1014707913437646E+07 

0.7189526088705342E+06 

0.5700357159533311E+06 

0.4852838429767232E+06 

0.4333653036170878E+06 

0.4003  407-5^78983  80E+06 

0.3792618208541028E■^06 

0.3663566904123807E-I-06 

0.3594433543536134E-t-06 

0.3572059996408999E-I-06 

0.3588425746327236E-t-06 

0.3638867262380186E-t-06 

0.3721195131488260E■^06 

0.3835335354846214E-t-06 

0.3983345162081381E-t-06 

0.4169781477992669E-t-06 

0.4402492810730193E-t-06 

0.4693969889513739E-t-06 

0.5063328348086139E-t-06 

0.5538348649144967E-h06 

0.6154173400618308E-t-06 

0.6936720638097440E-I-06 

0.7848578292114019E-t-06 


O  J  ,  Z,  -  J 

0  .  V  J  C  r.  -  ^  0 

0 . OCOCE-00 

O.OOOOE+CO 

0  .  OOOOE-^CO 

0 . COOOE+CO 

0  .  OOOOE-i-00 

0 . OOOCE-QO 

0  .  OOOCE-^CO 

0 . 0000E-*-G0 

0 . COOOE+00 

0  .  OOOOE-^09 

0 . OOOOE^CO 

0  .  OOOCE-f-00 

0  .  OOOOE-^00 

0  .  OOOOE-i-OO 

O.OOOOE-t-00 

O.OOOOE+OO 

0  .  OOOOE^-00 

O.OOOOE+00 

O.OOOOE-t-00 

O.OOOOE-hOO 

0 . 0000E-^00 

0 .  OOOOE-^00 

0  .  OOOOE-kOO 

0  .  OOOOE-hOO 

O.O0OOE■^OO 

O.OOOOE-fOO 

O.OOOOE-hOO 

0.00OOE-^0O 

O.OOOOE-t-OO 

0.0000E-^00 

O.OOOOE■^00 

O.OOOOE-t-OO 

O.OOOOE-t-OO 

O.OOOOE-t-OO 

O.OOOOE-^00 

O.OOOOE-t-OO 

O.OOOOE-t-OO 

O.OOOOE-kOO 

O.OOOOE-^00 

O.OOOOE-fOO 

O.OOOOE-kOO 

0.0000E-^00 

O.OOOOE-^00 

O.OOOOE-t-OO 

O.OOOOE-t-OO 

O.OOOOE-t-OO 

O.OOOOE-t-OO 

O.OOOOE-^00 

O.OOOOE-t-OO 

O.OOOOE-t-OO 

O.OOOOE-t-OO 

O.OOOOE-t-OO 

O.OOOOE-hOO 

O.OOOOE-t-OO 

O.OOOOE-t-00 

O.OOOOE-t-OO 

O.OOOOE-t-OO 

O.OOOOE-kOO 

O.OOOOE-t-OO 

O.OOOOE-t-00 

O.OOOOE-^00 
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0  .  ■5  53  375;607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0 . 6639752607619130E+02  - 
0.6639752607619130E+02  - 
0 . 6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0. 6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 


0 . 370000CCE+C0 
0 . 36C0000CS+00 
0.3500000CE+00 
0.340000CCE+00 
0.33000000E+00 
0.320000CCE+00 
0.3100000CE+00 
0.30000000E+00 
0.29000000E+00 
0.2800000CE+00 
0.27000000E+00 
0.2600000CE-t-00 
0.2500000CE+00 
0.24000000E+00 
0.23000000E+00 
0.22000000E+00 
0.2100000CE+00 
0.20000000E-HOO 
0.19000000E+00 
0.1800000CE+00 
0.17000000E+00 
0.16000000E+00 
0.15000000E+00 
0.14000000E+00 
0.13000000E+00 
0.12000000E+00 
O.llOOOOOOE+00 
O.lOOOOOOCE+00 
0.90000000E-01 
0.80000000E-01 
0.70000000E-01 
0.60000000E-01 
0.50000000E-01 
0.40000000E-01 
0.30000000E-01 
0.20000000E-01 
O.lOOOOOOOE-01 


3. 37:7333989994004E^C6 
0.9427755439230207E+C6 
0.1019815C82250440E4-07 
0.11^24278682^3669E■^07 
0. 177301492931178CE+07 
0.60907849346729C5E+07 
0.S699783983433943E+07 
0.2214935990570053E+07 
0. 159611339436071CE+07 
0.1525215336211312E+07 
0.1607614311235352E+07 
0.1793873295695752E+07 
0.2104183966a67282E+07 
0.2523179777238031E+07 
0.2940400301811221E+07 
0.3354138243385256E+07 
0.3502321774989490E+07 
0.3142732036521640E+07 
0.2690614618297547E+07 
0. 238564754443 7100E+07 
0.2212265301879473E+07 
0.2132126019325340E+07 
0.2114527588065502E+07 
0.2111049172969458E+07 
0.2072979041079385E+07 
0.2038270863269938E+07 
0.2130590848089927E+07 
0.3064759482032674E+07 
0.7261293817461066E+07 
0.5315399848594837E+07 
0. 233997 42790593 02E+ 07 
0. 1602618136046374E+07 
0.1422342783974448E+07 
0.3426647341610848E+07 
0.1883807543007345E+07 
0.1036305339370152E+07 
0.8958213026219722E+06 


0 . ccc ;e*co 

0 . C3C ;e-:3 
0 . CCC3E^G3 
O-w^jCE-^uO 
0 . OOOOE+00 
0 . OOOCE+00 
0 . OOOOE+00 
0 . OOOOE+OO 
0 . OOOOE+OO 
0 . OOOOE+OO 
0 .OOOOE+OO 
0 . OOOCE+00 
0. OOOOE+OO 
0. OOOOE+OO 
0 .OOOOE+OO 
0 .OOOOE+OO 
0 .OOOOE+OO 
0 . OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
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APPENDIX  C 


RHORES.IN  and  RES2.0UT 


b 

0 

(1 .25663706144D-6. 0 .DO) 

(5. DO, 0 .DO) 

(1 .25663706144D-6, 0 .DO) 
(8.85418782D-12, O.DO) 

0 

100 

(66.3975260761913D0, -.63D0) 
(66.3975260761913D0, -.63D0) 
1 


31 


30 


RHO  = 

(66.3975260761913, -0. 

630000000000000) 

ORDER  = 

0 

972  . 
ORDER  = 

2248651 

1 

972.2248651 

187.6500753 

1 S  7 . 5  5  C  0  7  5  3 

934. 
ORDER  = 

4283634 

2 

1906 . 653223 

195.4329582 

383.0830335 

984  . 
ORDER  = 

3801688 

3 

2891.033397 

185.8901814 

568.9732149 

915. 
ORDER  = 

8408388 

4 

3806.874236 

200.4964131 

769.4696230 

1013 
ORDER  = 

.711925 

5 

4820.586161 

131.9734826 

951 . 4431106 

884. 
ORDER  = 

6268189 

6 

5705.212979 

210.1082728 

1161.551383 

1062 
ORDER  = 

.136542 

7 

6767.349522 

176.3922360 

1337 .943619 

844. 
ORDER  = 

7351624 

a 

7612.084684 

224.9893403 

1562 . 932960 

1131 
ORDER  = 

.011708 

9 

8743.096392 

169.9882173 

1732.921177 

801. 
ORDER  = 

8185456 

10 

9544.914937 

245.6614577 

1978.582635 

1217 
ORDER  = 

.253981 

11 

10762 . 16892 

163.9662347 

2142.548869 

762  . 
ORDER  = 

8165370 

12 

11524.98546 

271.1924904 

2413.741360 

1305 
ORDER  = 

.358288 

13 

12830.34374 

159.8700718 

2573.611432 

735. 
ORDER  = 

4863637 

14 

13565.33011 

296.6414655 

2870.252897 

1358 
ORDER  = 

.412309 

15 

14924.24242 

159.6115646 

3029.864462 

728. 
ORDER  = 

4198871 

16 

15652.66230 

310.9098915 

3340.774353 

1327 
ORDER  = 

.425915 

17 

16980.08822 

165.7111766 

3506.485530 

752. 
ORDER  = 

0057504 

13 

17732.09397 

302.0255827 

3808.511113 

1197 
ORDER  = 

.332567 

19 

18929.42654 

181.9250964 

3990.436209 

820. 
ORDER  = 

1300882 

20 

19749.55662 

271.0676079 

4261.503817 

1017 
ORDER  = 

.080151 

21 

20766.63677 

214.2704989 

4475.774316 

949. 
ORDER  = 

4312538 

22 

21716.06803 

233.8585674 

4709.632883 

854. 
ORDER  = 

9080545 

23 

22570.97608 

270.90^201 

4980.540103 

1140 
ORDER  = 

.459325 

24 

23711.43541 

205.9048081 

5186.444911 

750. 
ORDER  = 

1737200 

25 

24461.60913 

351.0899253 

5537.534837 

1307 
ORDER  = 

.770698 

26 

25769.37983 

195.5022884 

5733.037125 

717. 
ORDER  = 

6688645 

27 

26487.04869 

406.8593350 

6139.896460 

1264 
ORDER  a 

.348619 

28 

27751.39731 

209.2200640 

6349.116524 

771. 
ORDER  = 

3806371 

29 

28522.77795 

371.4091562 

6720.525680 

1024 
ORDER  = 

.466302 

30 

29547.24425 

260.8257408 

6981.351421 

935. 
ORDER  = 

7814954 

31 

30483.02574 

295.8745499 

7277.225971 

32 


31 


806 . 50C3321 

31289 . 52663 

377 . 0388154 

7654.;c47;5 

ORDER  = 

1137.219276 

32 

32476.74595 

251.2729102 

7905  .537696 

'"'^707.3269736 

33 

33184.07293 

535.0537958 

3440.591492 

ORDER  - 

1257.312913 

34 

34441.38584 

257.4043703 

8697.995362 

ORDER  = 

741.1712328 

35 

35182.55707 

517.8036111 

9215.799473 

ORDER  = 

993.8564614 

36 

36176.41353 

339.7936620 

9555.593135 

ORDER  = 

931.8798625 

37 

37108.29340 

381.5471477 

9937 . 140283 

ORDER  = 

752.3541347 

38 

37860.64753 

572.7077589 

10509.84804 

ORDER  = 

1196.983720 

39 

39057.63125 

324.2430500 

10834.09109 

ORDER  = 

684.0315444 

40 

39741.66279 

824.7602214 

11658 . 85131 

ORDER  = 

1095.362600 

41 

40837.02539 

384.2368843 

12043 . 08820 

ORDER  = 

802.9128895 

42 

41639.93828 

613.8301144 

12656.91831 

ORDER  = 

783.3037057 

43 

42423.24199 

685.2079883 

13342 . 12630 

ORDER  = 

1084.661147 

44 

43507.90314 

455.9571781 

13798.08348 

ORDER  = 

656.6129334 

45 

44164.51607 

1334.284575 

15132.36805 

ORDER  = 

1065.854660 

46 

45230.37073 

526.9875837 

15659.35564 

ORDER  = 

752.4465220 

47 

45982.81725 

980.8526343 

16640.20827 

ORDER  = 

741.0281641 

43 

46723.84542 

1111.912906 

17752.12118 

ORDER  = 

1024.268083 

49 

47748.11350 

687.2612845 

18439.38246 

ORDER  = 

619.8824995 

50 

48367.99600 

2816.917682 

21256.30014 

ORDER  = 

933.2215920 

51 

49301.21759 

933.4278782 

22189.72802 

ORDER  = 

746.7748971 

52 

50047.99249 

1563.762956 

23753.49098 

ORDER  = 

633.1649123 

53 

50681.15740 

3453.033284 

27206.52426 

ORDER  = 

958.2016636 

54 

51639.35906 

1281.^03^856 

28487.56012 

ORDER  = 

582.4799094 

55 

52221.83897 

7948.997769 

36436.55789 

ORDER  = 

707.8216956 

56 

52929.66067 

3371.731383 

39808.28927 

ORDER  = 

760.9429371 

57 

53690.60360 

3126.318710 

42934.60798 

ORDER  = 

509.7821808 

58 

54200.38579 

39212540.21 

39255474.82 

ORDER  = 

735.8181638 

59 

54936.20395 

7248.565280 

39262723.33 

ORDER  = 

558.5152406 

60 

55494.71919 

9834.353822 

39272557.73 

ORDER  = 

462.4843988 

61 

55957.20359 

12355.77148 

39284913.51 

ORDER  = 

630.7404790 

62 

56587.94407 

99562.87663 

39384476.38 

ORDER  = 

398.3959575 

63 

56986.34002 

2346.981242 

39386823.36 

33 


32 


CRCER  =  64 

383.0784337 

CRDER  =  65 

57369.41846 

420.5474146 

ORDER  =  66 

57789.96587 

245.2046028 

ORDER  =  67 

58035.17048 

235.8824654 

ORDER  =  68 

58271.05294 

187.5687740 

ORDER  =  69 

58458.62172 

96.42096173 

ORDER  =  70 

58555.04268 

77.08895275 

ORDER  =  71 

58632.13163 

41.11082041 

ORDER  =  72 

58673.24245 

17.01825856 

ORDER  =  73 

58690.26071 

10.48916236 

ORDER  =  74 

58700.74987 

3.979771329 

ORDER  =  75 

58704.72964 

1.338944808 

ORDER  =  76 

58706.06859 

0.6696262256 

ORDER  =  77 

58706.73821 

0.2039898478 

ORDER  =  78 

58706.94220 

0.5706777265E-01 
ORDER  =  79 

58706.99927 

0.2417996204E-01 
ORDER  =  80 

58707.02345 

0.6470666331E-02 
ORDER  =  81 

58707.02992 

0.1500585757E-02 
ORDER  =  82 

58707.03142 

0.5436122764E-03 
ORDER  =  83 

58707.03197 

0.1382514275E-03 
ORDER  =  84 

58707.03210 

0.2642533004E-04 
ORDER  =  85 

58707.03213 

0.7976083188E-05 
ORDER  =  86 

58707.03214 

0.2077497211E-05 
ORDER  =  87 

58707.03214 

0.3344830613E-06 
ORDER  =  88 

58707.03214 

0.7956812163E-07 
ORDER  =  89 

58707.03214 

0.2183178614E-07 
ORDER  =  90 

58707.03214 

0.3255914900E-08 
ORDER  =  91 

58707.03214 

0.5763290314E-09 
ORDER  =  92 

58707.03214 

0.1515927923E-09 
ORDER  =  93 

58707.03214 

0.2542666128E-10 
ORDER  =  94 

58707.03214 

0.3342033362E-11 
ORDER  =  95 

58707.03214 

0.6909957503E-12 
ORDER  =  96 

58707.03214 

2331.335413 

393851i4.“’5 

1763 .133344 

39390917  .33 

470.7230732 

39391388.61 

441.9821828 

39391830 . 59 

197 . 3234042 

39392027.91 

88.53890809 

39392116 . 45 

60.67678287 

39392177 . 13 

21.31919876 

39392198 . 45 

10.13661601 

39392208.58 

4.949877259 

39392213.53 

1.429882639 

39392214.96 

0.6163998092 

39392215.58 

0.2336347882 

39392215.81 

0.5703041639E-01 

39392215.37 

0.2140896778E-01 

39392215.89 

0.6943875875E-02 

39392215.90 

0.1448718199E-02 

39392215.90 

0.4632894529E-03 

39392215.90 

0.1393918747E-03 

39392215.90 

0.2511985089E-04 

39392215.90 

0.6616063731E-05 

39392215.90 

0.1948834333E-05 

39392215.90 

0.3166628450E-06 

39392215.90 

0.6563098797E-07 

39392215.90 

0.1875309857E-07 

39392215.90 

0.3067539013E-08 

39392215.90 

0.4838736526E-09 

39392215.90 

0.1205394753E-09 

39392215.90 

0.2319829132E-10 

39392215.90 

0.2896975372E-11 

39392215.90 

0.5388991059E-12 

39392215.90 

0.1232142292E-12 

39392215.90 

34 


33 


0 . 1432734719E-12 

58707 . 03214 

0 . 1530812715E-13 

39392215 . 90 

ORDER  =  97 

0.1722497837E-13 

58707 . 03214 

0.1948379647E-14 

39392215 . 90 

ORDER  =  98 

0 .2411439458E-14 

58707 . 03214 

0.3921698957E-15 

39392215 . 90 

ORDER  =  99 

0.5155239694E-15 

58707 . 03214 

0 . 5821199916E-16 

39392215.90 

ORDER  =  100 

0.7977617019E-16 

58707.03214 

0.6797385709E-17 

39392215.90 

35 


34 


APPENDIX  D 


RHORES.IN  and  RES 1. OUT  (2nd  run) 


b 

0 

(1.256637061440-6,0. DO) 

(5.00,0.00) 

(1.256637061440-6,0.00) 

(8.854187820-12,0.00) 

58 

58 

(66.397526076191300, -.63478670700) 
(66.397526076191300, -.63478670800) 
1000 


36 


35 


Rho's  with  Functions  of  Rho's 


0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0. 6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0,6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0,6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  ~ 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 


0 . 6347867071019999E+00 
0. 6347867071030000E+00 
0.6347867071040000E+00 
0.6347867071050000E+00 
0.6347867071059999E+00 
0.6347867071069999E+00 
0.6347867071080000E+00 
0.6347867071090000E+00 
0.6347867071100000E+00 
0.6347867071109999E+00 
0.6347867071119999E+00 
0.6347867071130000E+00 
0.6347867071140000E+00 
0.6347867071150000E+00 
0.6347867071159999E+00 
0.6347867071169999E+00 
0.6347867071180000E+00 
0.6347867071190000E+00 
0.6347867071200000E+00 
0.6347867071209999E+00 
0.6347867071219999E+00 
0.63 47 8 67 07 12 3 000 OE+00 
0.6347867071240000E+00 
0.6347867071250000E+00 
0. 6347867071259999E+00 
0.6347867071269999E+00 
0.63 47867 07 12 80000E+00 
0.6347867071290000E+00 
0. 63 4786707 13 0000 OE+00 
0.63 47 867 07 13 0999 9E+00 
0.6347867071319999E+00 
0.6347867071330000E+00 
0.6347867071340000E+00 
0.6347867071350000E+00 
0.6347867071359999E+00 
0. 63478670713 69999E+00 
0.6347867071380000E+00 
0.6347867071390000E+00 
0.6347867071400000E+00 
0.6347867071409999E+00 
0.6347867071419999E+00 
0.6347867071430000E+00 
0.6347867.071440000E+00 
0.6347867071450000E+00 
0.6347867071459999E+00 
0.6347867071469999E+00 
0.6347867071480000E+00 
0.6347867071490000E+00 
0.6347867071500000E+00 
0.6347867071509999E+00 
0,6347867071519999E+00 
0. 63 4786707153 OOOOE+00 
0.6347867071540000E+00 
0.6347867071550000E+00 
0.6347867071559999E+00 
0.6347867071569999E+00 
0.6347867071580000E+00 
0.63 4786707 1590000E+00 
0.6347867071600000E+00 
0.6347867071609999E+00 
0.6347867071619999E+00 


0.99296 02 783900407 E+230.0CE+ CO 
0.1014410959638363E+240.0CE-00 
0. 103  63  472  602  3  53  30E+240 .OCE+00 
0.1059242052408832E+240.00E+00 
0. 1082871584885779E+240 . OOE+00 
0 . 1105445852  407965E+2  40. OOE+00 
0.1130 487245893498E+2 40. 00E+ 00 
0. 1156718772 688488E+2 40 .OCE+00 
0. 11836216572  08496E+24  0.00E+00 
0.1211302275900983E+240.0CE+00 
0.1240279174983865E+240.00E+00 
0.1270096260587070E+240.00E+00 
0. 1300324712146113E+2  40.00E+00 
0.1332135243075332E+240.COE+00 
0.13 655764324863 91E+2 40. OOE+00 
0. 13995611425 46825E+240.00E+00 
0. 1435776 62 0485633E+240.00E+00 
0.147288060 8 86 6953E+2 4 O.OOE+00 
0.1511678389003229E+240.00E+00 
0.15522 14 172 196980E+2 40. OOE+00 
0.1594141517130218E+240.00E+00 
0.163 8 082875 8 1803 4E+2 40. OOE+00 
0. 1682 98 82 8 4515 67 lE+2 40 . OOE+00 
0. 1730544418193773E+240 . OOE+00 
0.17 80372 5 83 4255 47E+240.00E+00 
0. 183 0 1 1053 42 6782 6E+240.00E+00 
0.18841922 1372 0950E+2 40. OOE+00 
0. 19403 66 83 025945 9E+240.00E+00 
0. 199917 0476567641E+240.00E+00 
0.2060998531770978E+240.00E+00 
0.2125112683449009E+240.00E+00 
0. 218790512 1564825E+240.00E+00 
0.2258244872771186E+240.00E+00 
0.2332320975378146E+240.00E+00 
0. 241057053 4841065E+2 40. OOE+00 
0. 24910658 80114460E+240.00E+00 
0.2577210069871970E+240.00E+00 
0. 266770053 1612090E+240.00E+00 
0.275 99 99962 8898 12E+240.00E+00 
0.2860212233708263E+240. OOE+00 
0.2964524343382726E+240.00E+00 
0.3076303921884608E+240. OOE+00 
'a.-'3193718259870463E+240 .  OOE+00 
0.3318334207454595E+240.00E+00 
0.3451108449759469E+240.00E+00 
0.3591539134708886E+240.00E+00 
0.3741651947534502E+240.00E+00 
0. 3 899233 65600 43 65E+2 40. OOE+00 
0.4068716235362964E+240.00E+00 
0.4247550501705975E+240.00E+00 
0.4436151269309913E+240.00E+00 
0.4641872996262146E+240.00E+00 
0.4861775406031703E+240.00E+00 
0.5096123701434384E+240.00E+00 
0.5351346910104587E+240.00E+00 
0.5624619966519761E+240.00E+00 
0.5911810151294268E+240.00E+00 
0. 620866769597693 8E+2 40. OOE+00 
0.6551448022409117E+240.00E+00 
0.6924004345741584E+240.00E+00 
0.7325167628672500E+240.00E+00 


37 
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0.6539752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0 . 6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0 . 6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0,6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
0.6639752607619130E+02  - 
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